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Abstract 

An eight-component formahsm is proposed for the relativistic two-fermion 
problem. In QED, it extends the applicabihty of the Dirac equation with 
hyperfine interaction to the positronium case. The use of exact relativistic 
two-body kinematics entails a CP-invariant spectrum which is symmetric in 
the total cms energy. It allows the extension of recent recoil corrections 
to the positronium case, and implies new recoil corrections to the fine and 
hyperfine structures and to the Bethe logarithm. 

PACS number: 03.65.Pm 



I. INTRODUCTION 

The relativistic two-body problem for two spin-1/2 particles is based on 16-component 
wave functions which transform as the direct product of two four-component Dirac spinors, 
7^(16) ^ ip^0^2- For unequal masses 1712 > mi, the equations are simplified by the elimination 
of the small components of particle 2 and by a subsequent power series expansion about the 
non-relativistic limit of this particle. One thus obtains an effective Dirac equation for particle 
1, with a hyperfine interaction that contains the Pauli matrices (T2 of particle 2. Such an 
equation has 4x2 = 8 components. It is very powerful for hydrogen and muonium 
More recently, a non-relativistic quantum-electrodynamics (NRQED) has been elaborated 
which allows one to eliminate the small components of both particles, which is particularly 
useful for the equal-mass case, mi = m2 as in positronium One thus arrives at 

an effective Schrodinger equation, in which the Pauli matrices cri and (T2 produce a four- 
component spin structure. However, the fact that the power series expansions in a are 
rapidly converging, does not prevent technical difficulties, presently at the order in the 
binding energies. These a^-terms have only been calculated by the above eight-component 
strategy, and only to first order in mi/ 1712. 

In this paper, a different eight-component equation is derived which does not eliminate 
small components and dispenses with nonrelativistic expansions. It exploits the fact that 
the chirality operator 75*71 commutes with the matrix 7^72 of the parity transformation; 
a corresponding separation of components does not exist for a single Dirac particle. We 
begin with the rederivation of an eight-component equation for two free spin-1/2 particles, 
in which the spin operator of particle 2 is removed in the cms (p^ = — = p) by 
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means of a new matrix c, given in ( ]20|) below. The spin dependence in the lab system is 
generated by a boost, to be discussed in section 5. The removal matrix c also mixes large 
and small components, via a Dirac matrix /3 = 7°. For total cms energy E, the result is 
an effective single-particle Dirac equation, for a free particle of reduced mass /i and reduced 
energy e: 

H = mim2/E, e = {E"^ - ml - Tnl)/2E, E = mi + m2 + Eb. (1) 

In section ^, the QED Born amplitude for cms scattering will be used to derive the interaction 
for this equation. The equation with interaction is for two leptons (no anomalous magnetic 
moments) 

{e-fiP- V)tlj = 75(0-1 - , = o-i X 0-2 . (2) 

Setting 75cri = ct = 7^7, the equation has the appearance of the usual Dirac equation with 
hyperfine interaction, particularly as ex. x (72/ E may be approximated by o; x (72/ {mi +7/12) 
to order a^. However, the complete operator —icx. x cr2V'p/E contains an antihermitian 
part which ensures exact relativistic two-body kinematics; for mi = m2, it produces the 
correct spin structure to order 0. A previous derivation from the 16-component Dirac- 
Breit equation produced a different hyperfine operator, which is equivalent to the new one 
only near threshold, —E^/E -C 1. The new hyperfine operator is left invariant by the c- 
transformation. (The previous derivation also had to assume a point Coulomb potential, 
V = —Za/r, a = e"^ , Ze = nuclear charge.) 

The relativistic on-shell two-body kinematics has long been well known. One has e^— /i^ = 
Ef — m\ = E2 — m\ = k'^ in the cms, where k"^ = —h? is negative for bound states, and 
the asymptotic form e"**"" of ip displays as multiple of a Bohr radius. The incorporation 
of relativistic two-body kinematics results in a spectrum which expresses E"^ in terms of m\ 
and m\. The combination Efi = mim2 is allowed in front of a square root (and possibly 
also with odd powers of Za, beginning with Z^a^). Our main nonperturbative result is 

E^ -ml-ml = 2mim2{l + Z^a^n*^)-^/^, n* = n - Pa , (3) 

where n* is an effective principal quantum number, and (3d is a quantum defect. Except 
for the details of f3d, (0) applies to any combinations of spins. It was first derived for two 
spinless particles [|l^], where also its empirical applicability to parapositronium was noticed. 
It was then extended to the case of one spin- 1/2 and one spinless particle The angular 



momentum defect, (3i = —61 = —{I' — I), has been discussed for two fermions of arbitrary 
magnetic moments to order Z'^a'^ It will be used in section 8 to derive a new "Barker- 
Glover" term in the fine structure. In section 4, a rather general formula for Z^a^ recoil 
terms will be derived which includes several new effects. In section 6, a vector potential 
is included in (Q) and evaluated in the dipole approximation. It leads to two additional 
quantum defects, 

n* = n + 6l-(3B-(3'6io, (4) 

where Pb is caused by the Bethe logarithm and (3' is an additional quantum defect in s- 
states. Parts of the Salpeter correction M are not included in M); their mass dependence 
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is examined in section 7. Vacuum polarization and nuclear charge distributions are also 
discussed in that section. For antiprotonic atoms, vacuum polarization must be included as 
a part of V for low-/-states. A method is proposed which extends the validity of to such 
states. However, also in cases where (0) must be solved numerically, the ii^^-dependence 
remains and adds new "recoil" corrections to the binding energies Ef, of (Q). 

It may be worth mentioning that the -E^-dependence of the spectrum is a very general 
consequence of the CPT-theorem. With the separate validity of C, P, and T in QED, one 
may also say that the -E^-dependence follows from C-invariance, but in the relativistic case 
the CP-transformation is slightly more convenient than C alone. Of course, the states which 
we calculate are CP-eigenstates only in the case of positronium. Muonium is transformed 
into antimuonium under CP. 



II. THE FREE EIGHT-COMPONENT EQUATIONS 

Let particles i{= 1, 2) satisfy the free Dirac equations, in units h = c = 1: 

{idi - OLiPi - mi/3i)'ilJi = , oci = 7^0-^, pi = 7° , 9° = d/dU . (5) 

The cTj are Pauli matrices, and 'jfPi + Pi'jf = 0. The direct product ip^^^^ = i^i ® "^2 
satisfies both equations © and thus also their sum, in which id^ + will be replaced by 
its eigenvalue , which is the total lab energy: 

(K° - 7i5piO-i - 72^p20-2 - mi/5i - m2/?2)V^^''^ = . (6) 

^(16) jg j^Q^ divided into two octets iI)lp and Xlp-, which have 7^ = 7!= 75 and 7f = — 7I = 
75, respectively: 

{K^ - 75PlCri - 75P2Cr2)V^LP = ("^lA + "^2/?2)Xlp , (7) 
{K^ - 75Picri + 75P2cr2)xLP = ("^lA + m2P2)ipLP ■ (8) 

The coupling between ipip and xlp arises because each Pi reverses the eigenvalue of •yf. In 
the chiral basis, 7^ and 7I are diagonal: 

7f-fi ?VA^f?il,^.^f^V^.P^fM,X.P^f^l, (9) 



'ii 



where the indices r and I (= righthanded, lefthanded) refer to the eigenvalues ±1 of and 
7|. In the parity basis, the Pi are diagonal, with eigenvalues ±1 for the large and small 
components g and /, respectively: 
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In this eight-component space, we also define a new matrix j3: 

/3 = /3i/32 = J , /375 + 75/3 = (12) 

(unit matrices are suppressed). Whereas [3i and 7^ do not commute, (3 does commute with 
7^7! , and both operators are diagonal in the representation (0) . As Lorentz transformations 
commute with 7^ and parity transformations commute with the decomposition of 
into if) LP and Xlp is invariant under the extended group of Lorentz (L) and parity (P) 
transformations. In the following xlp will be eliminated. We introduce a compact notation, 

P± = PiCTi ± P2cr2 , m± = m2±/5mi, (13) 

and observe P24'lp = Xlp, I32Xlp = i^LP in the basis ([TT|), such that we may effectively 
put = 1 and Pi = P in (0), (|) in this basis: 

{K'^ - -f5p+)iljLP = m+XLP , (K^ - 75P-)xlp = m+ijjLP ■ (14) 

Using the first equation for the elimination of Xlp, one obtains for the second one 

(ir° - 75p_)(m+)-i(ir° - 75P+)V^LP = m+^LP • (15) 

Multiplying this equation by m+ and using 



m_Lm_ = m 



m+75 = 75m_, (16) 



one arrives at the following equation 

ICo^pLP = , /Co = (ir° - 75l>- m_/m+)(irO - 75P+) - . (17) 

Equations (|1^ or (0) can be Lorentz transformed to the cms, where one has = 
E, = — = . The constants of (P!7| ) are combined into E"^ — m\ — m\ = 2Ee, with 

e defined by (P: 

[2_Ee - _E75(pcrm_/m+ + pAcr)] ^/^iP = 2mim2/3-?/'LP, cr = cri + 0-2 , Act = cri - cr2 , 

(18) 

where mim2 may also be written as according to ([^). If per were present without the 
factor m_/m+, one would use per + pAcr = 2pcri, and equation (p^Sf ) would be identical 
with equation (||) for V = 0, multiplied by 2E. 

It is in fact possible to remove m^/m^ from (|T^), by a transformation 

ipLP = cip , c~"^75 = 75C, cere = (Tm+/m^ , cAac = Acr . (19) 

Explicit forms of c are 

c = (m+m.)"^/^ 1712 + 1^1/5(1 + cricr2) = {Tn^'m_)~^^'^ {m_^. — 2miAs) , (20) 

where = (1 — cricr2)/4 is the projector on singlet spin states. To verify (p!9|), one notices 
AsCr = crAs = . In summary, the 16— component equation (^ is now transformed into a 
single free-Dirac equation (e — fiP — 75crip)?/'jree = 0, with no trace of the spin operators of 
particle 2. In the next section, it will be seen that the interaction for this equation contains 
no mass factors at all. 
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III. THE ONE-PHOTON EXCHANGE INTERACTION 



A connection between bound states and perturbative QED rests on the S'-matrix S = 
1 + iT and the Born series for the T-matrix, T = T^^^ + T^^) + . . . . When this series is 
summed by appropriate differential or integral equations, the bound states appear as poles 
of T. In detail, one takes plane waves, ipi = UiC'^'^' with 0j = kjrj — Eit for the initial states, 
0^ = k^Fj — E^t for the final states and extracts the resulting energy-momentum conserving 
5-function from the T-matrix elements, Sif = i{27r)'^6{E — E')6(ki + k2 — k'^ — k2)Tjj. In 
analogy, we put ip^p = ve^'^, Xlp = we*"^ with = 0i -|- 02- The first Born approximation 
to t/j of the matrix elements Tij is (with gi = — e , q2 = Ze , 



a) 



t-i) AirZa ,+ ,+ , , AnZa, ,+ , , ,+ , , , , , 

Tlf = ^— - (Xi(X2)UiU2 = ^^(1 - 0-10-2)^; + W'\l + (T,(T2)W] , (21) 

t t 

with t = q^"^ — , = — K^' , q = ki — k'^. Using (|14D for the elimination of f in one 
term and of w in the other, T-^^ is expressed in terms of an 8 x 8-matrix M: 

w = —(K^ - -i^k+)v , v'^ = w'Uk''' - 75A;'_)— , t/.'^ = w'^ Mv . (22) 

In the differential equation approach based on (^, the potential V is the Fourier transform of 
T^^\ apart from corrections from the hermitian part of T^^^ Unitarity S'^S = 1 implies 
rp{i)^ = T(i) and thus V = V\ i.e. a hermitian potential. However, the asymmetric form 
w'^^Mv implies a non-hermitian interaction JCj in our differential equation (/Co + /C/)?/' = 0. 
The simultaneous validity of (/Co + ^\)x = guarantees real eigenvalues, though. The 
ordinary Dirac equation fails to order Z'^a'^ in the states with I = f [l = orbital, / = total 
angular momentum), due to hyperfine mixing. 

We now restrict ourselves to the cms, ki + k2 = k'^ + ^ = and call ki = k , Pi = p = 
— iV. The total phase contains k(ri — r2) = kr, such that one has p2 = —p. Moreover, 

= = E. To order Z^a^, one also has g° = , t = — q^: 

- q^M/AnZa = {E - -i^k'_)m-\l - a^a^) + m-\l + (t^ct2){E - 75^;+) . (23) 

With k'_ = k'cr, one has k'_{l — a 1(72) = 4k'crAs = 0, as 1 — (Ti(T2 vanishes for triplet states, 
while cr annihilates the singlet state. Consequently, 

- q^M/inZa = m-^[2E - 75(1 + <Ticr2)kAcr] = 2m;^(E - i75k<T'') . (24) 

Combining now m^^ with w'~^, the Fourier transform of m^M will be called JCf. 

JCi = EV - t-f^Vpcr'' , V = -Za/r. (25) 



This is to be added to /Cq (|T7|) and used in (0). To arrive at the form (^ of the differential 
equation, /C/ must be replaced by c~^/C/c according to (plQl). Fortunately, one finds 

ccr^c = cr^ c-^/C/C = /C/ , (26) 

such that the same operator appears in fact in (^). To order Z^a^ in E, /C/ is equivalent 
to the operator )Cjb = EV — '^^{icr^ — 0-2) [y,p] of the Dirac-Breit approach [§,0, the 
transformation being 
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(27) 



after the transformation c (notice c^^Kib c ^ Kib)- 

For the orders Z^a^ and Z^a^, the condition g'' = should be replaced by the current 
conservation conditions, g^Jf = {i = 1,2). This produces the spin operator cr;j'-cr^ = 
fiC2 — {<^iQ){(^2<i)/q'^, which is equivalent to using the Coulomb gauge. However, we have 
not yet performed this calculation. The new Z^a^-terms of (^) below result from the use 
of relativistic kinematics in (|^). 



IV. SOLVING THE NEW EQUATION NONPERTURBATIVELY 

The total angular momentum F = L + (cri + (T2) /2 = Ji + cr2/2 is conserved. All eight 
components of have the same eigenvalues /(/ + 1) = F"^ of F^ and rrif of F^; factors 1/r 
and i/r are extracted from ipg and ipf as usual: 

i^9k = Xk^^'^Ugkir) , ^fk = xi'^'^-Ufk{.r) , (28) 



[e — ji — V)ug = inuf , (e + /i — V)uf 



-ITTU 



a ' 



TT 



[(Ti-i(T'^V/E)r^/r . (29) 



The index k assumes the four values /c = s = (/, 0) (singlet). A; = 1 = (/, 1) (triplet with 
orbital angular momentum / = /) and /c = + , k = — (triplets with / 7^ / and eigenvalues 
=F1 of <Jir<72r = (Tif(T2r) M. Wc arrange the radial wave functions as follows: 
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(30) 



with d± = dr±l/r, /± = 1±2V/E. The states k = ± with / = /±1 have ji =32 = /±l/2 
to order a^. An analytic solution of the equations is obtained when the hyperfine operator 
is replaced by an equivalent r~^-operator. At this level of precision, the equation becomes 
equivalent to the one derived from the Dirac-Breit equation [^,|^, where the replacement 
was achieved by the substitution r = r' — Zaa/2fi (see section 8 for a). 



0. 



(31) 



LF' = L?' — Z'^a^il + a) comprises all r ^-operators. Its eigenvalues will be denoted by 
/'(/' + 1), to profit from the analogy with the Schrodinger equation. Defining moreover 



-2ikr = 2Kr, az = Za, az^/i^ 



n 



(32) 



and u = e ^^'^z^' F{z), the equation becomes as usual an equation for the confluent hyperge- 
ometric function F, 



[/' + 1 - zdl - {21' + 2- z)d,]F = n*F . 



(33) 
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For unbound states, t] = —in* is called the Sommerfeld parameter. For bound states, 
n* — I' — 1 must be a non-negative integer n^. With /' near an integer /, one defines 



61 = 1' -I 



n 



rir + 1' + 1 = n + 61 , 



(34) 



where n = + I + 1 is the principal quantum number and 61 is always negative. All four 
eigenvalues A of — + q;| are given by 



A"[A' - 2A - 4(F" - a'ze/E) + 2ai.] - 2a'z\ 



;i -4e/E) = 0. 



(35) 



It contains the Dirac eigenvalues with recoil-corrected hyperfine structure and hyperfine 
mixing near both static limits 0. For equal masses, the equation is reliable only to order 
a^, but for / = / we nevertheless quote the result for 61 to order a"^. With l—ie/E = Q;^/4n^, 
one solution of (|35| ) is A = — a^/8n^. It belongs to parapositronium, as we shall see. The 



value of /' follows from I' + \ 



F2 + 1/4 + A - a2 



-61 



a 



para 



+ 



a 



a 



21 + 1 (2/ + 1)3 8/22(2/ + 1)' 



(36) 



For orthopositronium with / = /, one factor A is divided off and A^ can be neglected to 
order a^: 



— ^hrthoi^ — f) 



a 



21 + 1 



1 



a 



2F2 2F6 



a 



1/2F2 



2/ + 1 



a 



8n2(2/ + l) 



(37) 



We now come to our main point, namely the calculation of E"^. From ( P^ and the definition 
(0) of /i and e, one finds (0). Expansion of the square root leads to a more practical series, 
which to order a| is 



E'^-m^ 

17111712 



n 



*2 



Sal 
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5 ai 
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Insertion of n* = n + 61 yields, to order a|. 



E' 



m 



mim2 



l + ai 



n . 



261 



a. 



a. 



m 



nil + ^2 • 



(38) 



— + 

An n 



a^n 



a. 



(39) 



This expression is still quite compact, in view of the fact that 61 contains both fine and 
hyperfine interactions. The result for E is, again to order a| 

E = [m^ — mim2a|(l + a\b/n)/n^Y^'^ = m — jjnrCt'zi^ + a\b/n)/2n^ 

- AiL«l(l + alh/nf/9>n^ni - ^3^a|/16nW, (40) 

where = inim2/m is the non-relativistic reduced mass, and the expression (|39| ) for b 
remains to be inserted. The third term of (|40|) , with the approximation 6 = 0, is known as 
the Bechert-Meixner recoil correction To order a|, b may be approximated by —261/a'^ — 
3/4n. For comparison with the literature |1T3|J^, we split the reduced mass finr in a rather 
unusual way. 
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fJ-nr = fJ'nril + fJ'nr/m) - filjm ^ ml{l - m\/m^) - filJm. (41) 

In the order a|, we may then combine all contributions from (40) with that of an additional 
operator ~ L^/r^ [|13|, which makes the expression complete for / > 0: 

The last contribution has been symmetrized in the masses, and the expectation value 
refers to the solution of the Schrodinger equation with reduced mass. This operator has been 
calculated from two-photon exchange, but the order of exchange may be gauge- dependent. 
Insertion of the Dirac quantum defect 

6ln/o.l = (7 - J - ^ -[1 + «l/(2j + l)1/(2j + 1) (43) 

reproduces the known result [^. Equation (^21) generalizes this result to arbitrary masses 
and hyperfine interactions (section 8). 

However, we would like to advocate the direct use of formula (|39|) for {E'^ — m^)/mim2, 
because it is more compact for calculations and less mass-dependent for measurements. 



V. BOOSTS AND COORDINATE TRANSFORMATIONS 

Lorentz transformations can be constructed directly for ip and ipLP, but it is more con- 
venient to use the known transformations of the single free-particle spinors ipi and iIj2- We 
only need the boosts from the cms to the lab system (/), where the system has a total 
four- momentum K^^: 

= AV'^cm, = (7 + llKaif'^ = (27 + 2)-^'\^ + 1 + ^iKcn) , (44) 

K = K/E, 7 = K7E= (l+i^2)i/2_ (45) 

We also take the 2;-axis along K, Kcji = Kaiz. Suppressing the index ^p, the eight- 
component lab spinor is 

ijl = A^err^, A = (7 + ^^ka2,f'\l + 75^^^^!.)'/' = (1 + \K^oI + i^a,775)^/' (46) 

A = 7 + \k^,a, - \k^Aal/{2^ + 2) (47) 

(in checking ( ]i7| ) by squaring of (^), use (Jz^az = 0, (Acr^)^ = 2(1 — crizcr22), (Acr^)^ = 
4(AcT2)^ and = -^'^(7^ — 1)). The boost A for x follows from A by replacing 750-2 by 

-15(^2- 

Xi = Axcm, A = 7 + Ik^^^az - \k^al/{2^ + 2). (48) 
The inverse boost has K replaced by —K, which is equivalent to a sign change of 75: 
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Af3 = pA-\ A(3 = pA-\ (49) 
Insertion of (19), ipcm = cip gives in the lab system ipi = Acip, and for the c-transformed ipf. 

^ic = c-^^i = c-^Acifj = A,ij, A, = c-^Ac, (50) 

^ = 7 + l^,Ka,m+/m_ - \k''^al/{2^ + 2). (51) 

The corresponding boost for x, on the other hand, has a factor m+ extracted: 

A^ = m-^Am+ = 7 + l-f^KAa^m+fm^ - \k^al/{2-i + 2). (52) 

The desired boosts for ifj and x are A^. and Ac, respectively. They are needed for the 
construction of Dirac-Breit equations in the presence of external potentials. 

In a covariant treatment, the interaction between two particles at distance ri = vi — r2 
depends also on a time difference = ti — ^2, such that = {x^, ri) is a four- vector. 
A second independent four-vector X'^ is defined such that = idx becomes the total 
four-momentum pi + P2, which is conserved: 

= - x^, = Eix'^ + Eix'i, pt = p^ + EiPf", pi^ = -pf" + i'aP'^, (53) 



^1 + ^2 = 1, P''=pt+pl P^ij = K^i^. (54) 

The as yet open value of Ei — E2 is chosen such that p° = id/dx^ vanishes at asymptotic 
distances in the cms where particles 1 and 2 are on their mass shells, Ef — k"^ = ml, E'^ — k'^ = 
i.e. E2 — Ef = m\ — m\ = m^m^. Using in addition Ei + E2 = E, one obtains 

El = (P2 _ m+m_)/2E, E2 = (^^ + m+m„)/2E. (55) 

Extracting now from ( |53D 

p>^ = E2pt - Eip'i (56) 

and inserting the asymptotic values E^ of p° in the cms, one finds that p^ vanishes here for 

Ei = EJE = i(l =F m+m_), m± = m±/E. (57) 

In the interaction region, p^ does not vanish, but in the context of a single integral equation 
(Bethe-Salpeter equation) it cannot be treated as a dynamical variable. Instead, integrals 
involving p^ are treated as perturbations on 5 (p°) -integrals (section 7). 
The space coordinate transformation of (|53D is 

ri = R + E2ri, r2 = R-Eiri, (58) 

where zi is Lorentz contracted, zi = 2/7. The nonrelativistic approximation yields the 
familiar E^ = mi/m. 
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VI. VECTOR POTENTIAL AND BETHE LOGARITHM 



In the presence of a four-potential A'^, and P2 are replaced by 

< = P^ + eA'^(xi), n!^=pi^-ZeA^^ix2). (59) 



The coefficients Ei and E2 of the transformation (|53D to vr^ and 11^ remain unchanged in a 
perturbative treatment of A^. In the dipole approximation A{xi) = ^(0:2), (|§) leads to 

77 = ^27ri - ^1772 =p + rcA, (60) 

where re = r ■ e is the dipole radiation charge [g| , 

r = i(Z + 1) - i(Z - l)m+m_ = 1 + i(Z - 1)(1 - m+m_). (61) 

In the following, we use the Coulomb gauge in the cms and replace p by tt in (0) . To exhibit 
its i?^-dependence, we also multiply the equation by E: 

{Ee -EV - mim2/3)V = 75(0-1 - «<T^1^/E)E7r^. (62) 

The substitution 

r = Ep, p = Pp/E, 7r = 7r,/E (63) 

brings (^) into a form which contains only even powers of E. For the moment, we assume 
V = —Za/r, EV = —Za/p = V{p). Deviations will be discussed in section 7. For the 
Lambshift calculation, one neglects the hyperfine operator and obtains an explicit eigenvalue 
equation for Ee = {E"^ — m\ — m\)/2, 

hip = Eeip, h = mim2(3 + V{p) + aiTp. (64) 

Its Coulomb Greens function G satisfies the usual equation [Q], taken in the new variables 
P, P'- 

[Vj + {Eef - mlml + 2EeZa/p + {Za + iap)Z'^a^ / p^]G{p, p', Ee) = 5{p - p'). (65) 

Evidently, it is also independent of the signs of nii and 1712. The same remark applies to 
r (|6l|) , but that expression uses the low-energy dipole approximation. Apart from that, 
the Bethe logarithm can only depend on E"^, m\ and m\. Its proportionality to r^ has 
been noted previously 0,^. Inspection of the formulas collected in [§] reveals another 
small s-state correction, which is also proportional to r^. Moreover, both corrections are 
proportional to mim2 and therefore pushed under the square root in (^), where they appear 
as quantum defects: 

PB = ^r'\nko{n,l), P' = -^r' - In a') . (66) 
After the extraction of Pb and j3', there remains a somewhat reduced Salpeter shift. 
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(67) 



-25, 



/o 



2a " 1 1 1 

In h>-H 

n ^ i 2 2n 



/(/ + 1)(2/ + 1)' 



(68) 



with a'-^ = — 21n(2a;) — 2. The denominator mj^m_ in (|67| ) will be discussed in the next 
section. The factor 1/m in front disappears in the expression for E"^, as E"^ m? + 2mEi, 
according to (0). 



VII. BETHE-SALPETER EQUATION, VACUUM POLARIZATION AND FORM 

FACTOR 

Continuing with the T-matrix of section 3, the two-photon exchange part T^^^ may be 
used to derive corrections to the main interaction, including the Salpeter shift (|U7|). However, 
here we consider instead the more commonly used Bethe-Salpeter equation. It is formulated 
in the cms, where ( ^3|) reduces to 

p^, = p'^ + E^g^o, p^ = -p" + E^g^o, (69) 

and p^ is an integration variable. Suppressing the integration, the BS-equation applies to 
the two-fermion Greens function G = G^^^^ 

G-dS)^ 5(16)^^(16)^(16)^(16)^ 

where S^^^^ is the product of two free fermion propagators in momentum space 

5(16) = (^^ _ ^^)-i ^ (^^ _ ^2)-i = (^^ + ^ (^^ + m2)/D,D2, (71) 

D^=pl-ml + ie = p^ + 2Eip^ + k"^ + ie, D2 = pi - ml + is = p^ - 2^2]?° + k'^ + ie, 

(72) 

and the kernel K^^^^ is the sum of all irreducible Feynman diagrams. Our proposed eight- 
component formalism simplifies 5* at the expense of K. Taking K from the matrix elements 
X^Kip analogous to the 8x8 Matrix M in (^2]), one has 

5(16)^(16) ^ 5(8)^(8)^ 5(8) ^ ^ VLpwlp/D^D2. (73) 

In a matrix notation where = (3 and 75 are identical for particles 1 and 2, the spin 
summation is now 

^vlpwIp = m2-Ei + miPE2 - m_p^ - 75p(m2<Ti - mi(3cT2). (74) 

It is only a first-order polynomial in pi, p2. Consequently, the numerator of S*^^-* is linear in 
the integration variables, whereas S**^^^-* is quadratic. One may in fact go one step further 
and use irreducible representations of the Lorentz group, in which case one arrives at S*^^-* = 
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l/DiD2. This amounts to the ehmination of the dotted spinor components from ijj |T^; it is 



widely used in QCD calculations ^7j. Four bound states, however, it has the disadvantage 
of suppressing explicit parity invariance. 

The last bracket in (ff^) is (m_cr + m+Acr)/2 in the notation (|l^), and m_cr may be 



converted into m+cr by the transformation (|20| ) 

Vlp = cv, Wlp = c'^w, ^ fw"'^ = (e — 75crip + /i/? + p°m_/m+)m+. (75) 

To handle the p^-integration, one may use the formula 1 / {x—ie) = P/ x+i7r6{x) (P=principal 
value) : 

2E/D1D2 = 2iix6{p^)/{p' - e) - l//^i(/ + le) + 1/D2(/ - le), (76) 

and treat all p^-integrals except the first one as perturbations of the Greens function. In 
these integrals, the spin summation (^) reduces to J2vwym+ = e — ctp + which is 
simply the expression for a single-particle spin summation. In the remaining integrals, the 
complete two-fermion propagator has the more explicit numerator 

2Ej2'^w^ /^+ = E"^ - ijnl + ml- 2mim2/5)(l - Ep^/m+m^) - 275^^0-1. (77) 

After integration, this p'^-dependence produces factors {ml + m^) /nij^m^ and 
mim2/3/m+m_. Such factors may also arise from the denominators I/D1D2 of the 
crossed graph, which contain the combination p^{Ei — E2) = Ep^m^rh_. To check the 
mass dependence in the 5;o-piece of ^E'^^^ (|67D, notice rn^lnlrn'^ /m2) — m^ln(m^/m^) = 
m^m_ \\i{'mim2 / m?) + {m\ + \ji{m? / 77111112) . 

We conclude with a discussion of the potential EVij-) V{p) in the presence of vacuum 
polarization and nuclear charge distribution. For a heavy particle 1 (yU~ or p), electronic 
vacuum polarization is so large in low-/-states that it must be added to the Coulomb potential 
Vc, in the form of the Uehling potential: 

V = Vc + Vu, Vu = -^ dX'e-^^S{X'), (78) 

Jvrr JO 

with = (1 - 4m2/A2)V2(i + 27nl/X'^)/X^Q{A77il) in the electronic part of Vu- At fixed A, 
the r-dependence is easily converted to a p-dependence, {E/r)e~^^' = p~'^e~^^P: 

Vu{p) = -^ r dX\-'^^S{X'). (79) 
Snp Jo 

The substitution XE = X' shows that the integral depends only on E"^. There is also a 
second-order recoil correction V^^J which is spin-independent [0. For antiprotonic atoms. 



the radius of the vacuum polarization cloud is much larger than the Bohr radius {Zafi)~^: 
below a critical value Ic of the orbital angular momentum /, vacuum polarization exceeds 
all relativistic effects (protonium pp and p^He have Ic = 3 and 7, respectively), for all 



values of n [IS]. For these orbitals, V is best constructed numerically from (ffOD, particularly 
for the calculation of annihilation which is localized at very small p. However, for the 
calculation of the fine and hyperfine structure of these inner orbitals, the form (^) suggests 
the introduction of a running electric coupling constant ae'- 
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V = Ve + 5Vu, Ve = -Zajp, 5Vu = Vu + Z{ae-a)/p, {5Vu) = Q. (80) 

In this manner, the nonperturbative result (3) remains apphcable, with > a. 

Recoil corrections to an extended nuclear charge distribution Pc{t) are particularly large 
for muonic atoms. 

V{p) = -Za J dh\p - r/^|-Vc(r), Pc{r) = J d'^qe'^l'' F{q'). (81) 
Again, V{p) is i?^- dependent. 



VIII. ANGULAR MOMENTUM DEFECTS AND BARKER-GLOVER TERM 



To order a| = Z'^a'^, 61 is easily found for particles of arbitrary masses rrii and (^-factors 



gi = 2(1 + Ki). 61 is conveniently expressed as follows: 



6l/al 



-(l + a)/(2Z + l). 



With the abbreviations mj = rrii/E, i = e/E and the combinations 



e + 2Kim2, 



TTlr, 



= 1 + 2/«im2 + 2/t2mi + 2e(l - 5'i5'2/4), 



one finds 



a(/ = / ± 1) = ±cV(2j + 1) ± gig2i/{^f 
(where a sign error of H] for / = / + 1 is corrected), and 



a{l = f) = ± 



4F2 



(2/+i: 



„2 2 

16F2£i^ T 1 



2c?c 



,2^2 
2 



2j + l c2(2/ + l) 



1 + 



(82) 



(83) 
(84) 



(85) 



AF^cjcj 
d{2f + 1)2 

(86) 



In (^51), / = / ± 1 implies j = ji = / ± | = / =F |- For (?! = (72 = 2, c = 1 leads exactly to the 
Dirac fine structure component ( |13| ) of 5/, for all values of mi and m2, including mi = m2 
as in positronium. The case / = / is more complicated, because the two values j = I + ^ 
and J = / — I are mixed by the hyperfine interaction. We therefore define j as that value in 
the integer 2j + 1 which appears in the main term c^/ (2j + 1) of a{l = f) after expansion of 
the square root. It is frequently said that the fine structure contains the "Barker- Glover" 
term (for / > 0) R,|], 



E 



BG 



«l/"'(^-j)/[(2j + l)(2/ + l)mV]. 



(87) 



This is true for a spinless nucleus, except that the electron's anomalous magnetic moment 
frequently reverses the effect |jlT|. On the other hand, the weighted average over the hyperfine 
structure produces a more complicated expression. 



(/ - j)[2c' + g,g,e/2 - 2(cfc^/c^)(l + AL'ci4/c\2l + l)2)]/(2j + 1 



(88) 



13 



This is the final result, and as a rule it is much larger than the Barker- Glover term. For 
Qi = 92 = 2, however, c\c\ ~ /i/m e to this order in a, and the main terms cancel in (p8|). 
There remains a small rest, 

= alfXnr5l/{2l + l)n^ = -/i„^a|a/n3 = 4/i^^a|(/ - j)LV(2/ + lf{2j + l)mV, (89) 

which replaces Ebg in the case of muonium. 
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